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We study the Einstein gravity and dust system in three spacetime dimensions as an example of a 
non-perturbative quantum gravity model with local degrees of freedom. We derive the Hamiltonian 
theory in the dust time gauge and show that it has a rich class of exact solutions. These include 
the Banados-Teitelboim-Zanelli black hole, static solutions with naked singularities and travelling 
wave solutions with dynamical horizons. We give a complete quantization of the wave sector of the 
theory, including a dehnition of a self-adjoint spacetime metric operator. This operator is used to 
demonstrate the quantization of dehcit angle and the fluctuation of dynamical horizons. 
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I. INTRODUCTION 

The difficulty in formulating a four dimensional theory of quantum gravity has led to the study of many simpler 
models. These include symmetry reductions of 4D general relativity [1] and dimensional reduction to 3D gravity [2-7]. 
There is a large volume of literature on the latter, which includes pure gravity with point defects and/or topological 
degrees of freedom [8-11], topologically massive gravity [12], and higher spin gravity. While some of these (lower 
dimensional) simplifications have allowed for covariant quantization, there is relatively little work on the canonical 
quantization of any gravity-matter model. 

Our purpose in this paper is to develop a 3D model of gravity with matter which has the potential for complete 
quantization. This would reveal insights into quantum gravity in a setting with local degrees of freedom. The 
pressureless dust matter we use is perhaps the simplest such model, but it is sufficiently non-trivial in that there is 
a rich class of classical solutions, including ones with dynamical horizons. Such solutions are of much interest at the 
quantum level: questions such as what is the quantum analog of a classical dynamical horizon remain unanswered, 
and are key to understanding what is a “quantum black hole” or a “quantum trapped surface” [13]. 

With these issues in mind we begin by formulating a canonical theory of 3D gravity coupled to pressureless dust. 
This is a special case of the Brown-Kuchar [14] model which is designed to provide a dynamical matter reference 
system for gravity in 3-1-1 dimensions. It was used to give a physical Hilbert space setting for loop quantum gravity 
in the dust time gauge [15-17], and added as an additional world sheet field in the bosonic string to yield a curious 
extension of that theory [18]. 

We will see that in the 2-1-1 model, the dust time gauge gives a physical Hamiltonian that describes the dynamics of 
one local geometry degree of freedom; this remains in the circularly symmetric setting we consider in detail. The model 
also provides an example of the transfer of a matter degree of freedom to a geometric one; this may provide a useful 
viewpoint for quantum gravity in a more general setting, distinct from the strict conventional separation of matter 
and geometry degrees of freedom. In Section H we develop the circular-symmetry-reduced theory and in Sec. HI we 
give the gauge fixed theory. In Section IV we give several types of classical solutions, followed by the construction of 
a quantum theory of the system in Section V, with focus on the travelling wave solutions. The concluding section is 
a summary of our results and a discussion of further questions. 


II. HAMILTONIAN THEORY 

In units where 8G = c = 1 the action for gravity and dust is a sum of the two components 

S = Sg + Sd. 


( 1 ) 
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Let us consider this action defined on a three-dimensional manifold with topology E x M. The gravitational part of 
the action is 

Sg = ^ J (2) 

where is the scalar curvature of spacetime and A is the cosmological constant. The dust action is 

Sd = -^ J , (3) 

where m{x) is a function of the spacetime coordinates. 

To derive the Hamiltonian formulation we use the ADM (Arnowitt-Deser-Misner) parametrization of the line element 

ds^ = -N^dt^ + qab{dx^ + N‘^dt){dx^ + N^dt), (4) 


where qab is the space metric, N is the lapse function and N°' is the shift vector. With this the gravitational part of 
the action takes the well-known canonical form (see eg. [6]) 


Sg = ^J dx^ {n-^qab - NU^ - N-C^) , 


(5) 


where 7f“^ is the (density weight one) momentum conjugate to qab- and 7V“ appear as the Lagrange multipliers 
corresponding respectively to the Hamiltonian and diffeomorphism constraints 


Cf = - 2 Vaif,“, 


where q = det qab, '^ = ^ad ^'^^R is the Ricci scalar of the spatial hypersurface. 
The canonical dust action is obtained starting with the momentum 

SSd 


Prh := 


which leads to 


Sd = 

= 

Cf = 


1 / 


S(j) 

[ dx^ 

Pi 


N 


-N‘^da 


2 \ rriy/q 


- Nn^ - , 

my/q{q°‘^da(l)db4‘ + 1 ) ) , 


P4,da(t>- 

The variation of m in the dust action gives the equation of motion 

Pa 


m = ±- 


Using this we eliminate m from the Hamiltonian by writing the dust part of the scalar constraint as 

= ±P,j,\/q°-^da4>db(l) + 1 . 


( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 


The sign ambiguity will be determined below when we fix a time gauge. 


A. Imposing circular symmetry 


Let us now impose circular symmetry. A parametrization of the ADM phase space variables {qab,^°'^) for circular 
symmetry may be prescribed by using the flat 2D metric Cabdx'^dx^ = dr^ -I- r^dO^ and the radial vector s“ = [1,0] 
and Sa = s^Cab = [1,0]. In these coordinates 


<lab{t,r) 


Vi^{t,r))SaSb P 


J .2 


{eab{r) 


SaSb) 


^2^ „ 6 

2D(t,r) 


r^Ppit,r) 

2p{t,r) 


(e“*'(r) - s“s'’) . 


(14) 

(15) 
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With these definitions the symplectic term in the gravitational action is 

^ J dr d9 dt TT°-^qab = J dr dt (^Ppp + ■ 

The ADM metric becomes 


(16) 


ds'^ = - (iV2 _ df + 2n^N^dr dt + n^dr^ + p^de"^, 

and the Ricci scalar on the slice is 


vip \ n ^ 

Adding the gravitational and dust parts, with ^/q = |Dp|, gives the symmetry-reduced action 


S = J dr dt (Ppp + PnSi + P^<j) -Nn- N^Cr) , 

^ = sgn(Dp) 1^2 (^0 _ +2A|Dp|±P^^^^) + 1, 

Cr = p' Pp — ^Pq + P4>4‘' ■! 


(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


where we have used ‘primes’ to denote derivatives with respect to the radial coordinate. As one would expect, the 
angular component of the diffeomorphism constraint is identically zero {Crf, = 0) so that only radial diffeomorphisms 
play a role in the symmetry-reduced theory. 

The Poisson algebra of the constraints is first class 

{H{N),H{M)} = Cr{NM' - MN’), (22) 

{Cr{N^),Cr{M^)} = Cr {N^M^y - M^N^Y), (23) 

{H{N),Cr{N^)} = -H{N'N^), (24) 


being the reduced version of the Dirac/ADM algebra. At this point, with gauge freedom remaining, there are three 
pairs of conjugate variables parameterizing the six-dimensional phase space. The physical theory obtained by a Dirac 
gauge reduction, which fixes the constraints and removes the gauge ambiguity, will leave only one pair of conjugate 
variables in the two-dimensional physical phase space. In the following we consider the case of non-compact spatial 
slices with full gauge fixing and appropriate boundary terms to obtain a well-defined variational principle for the 
canonical 2-1-1 action. 


III. GAUGE FIXING AND PHYSICAL HAMILTONIAN 


In this section our goal is to obtain the Hamiltonian theory of the local physical degrees of freedom by fixing gauges 
and solving the constraints. 

We first fix the radial coordinate gauge by imposing Xp •= p — r 0. This is a standard choice in spherical 
symmetry, and is second class with the diffeomorphsim constraint 

{Xp,CriN’-)} = NP (25) 

Keeping this constraint preserved dynamically gives a relation between the lapse and shift functions, 

iV’'-sgn(D)^=0. (26) 

Solving the diffeomorphism constraint and imposing this gauge condition removes p and Pp from the system. We have 

Pp = npl, - P^c^', (27) 

which leads to the the partially gauge-fixed action 

S = J dr dt (^Pnfl + P^<y - Nn^ , (28) 

n = isgn(D)Po - nPl,) + 2A|D|r + ± (29) 
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In writing the Hamiltonian we assume that sgn(H) is fixed since we must have |H(t,r)| > 0 for the metric be non¬ 
degenerate. 

We now choose the dust time gauge by adding the constraint x<l> '■= (p — t 0, a. condition which is second class 
with the Hamiltonian constraint: 


{X^,H{N)} = ±N. (30) 

Requiring that this gauge is dynamically preserved leads to 

7V = ±1. (31) 

Recalling now the definition of the momentum (8), we see that the signs of m and N are linked in this gauge by 
rriy/q = NP^j,. Therefore choosing iV = 1, which generates dynamics forward in time, fixes [19] the sign ambiguity 
arising from (12): = +my/q and the solution of the Hamiltonian constraint gives the physical Hamiltonian density 

-Hphys = P4.= - 2A|H|r - . (32) 

The shift function is also fixed via (26): 

N=+l ^ = sgn(H)^. (33) 


A. Reduced action 


The reduced action is obtained by substituting the gauge fixing conditions and the solutions of the constraint (32) 
into the starting action (19). This gives 


where 


S = 


dt f dr (^PqCi - "Hphys) 

J 2 



^phys 


2A|H|r- \n\ 



(34) 


(35) 


The boundary term arises from the total derivative present in (32), and comes ultimately from the Ricci scalar density 

^(2)i? = (2/|H|)'. _ 

If St is asymptotically flat, as will be the case for some but not all solutions to the equations of motion, this term 
evaluated at a fixed radius 


2 


|H| 


(36) 


determines the energy within a disc of radius tq [20], and as we shall see below, this term also gives a measure of deficit 
angle at the origin in the limit tq —>■ 0. This is because in 3D gravity a conical defect represents a point source of 
energy at the origin; the relationship between deficit angle a and the energy M of the point source (in units SG =1) 
is M = a/TT [2]. 

Many interesting solutions in 3D gravity are singular at the origin [5], and in order to allow for these solutions 
we excise the origin r = 0. This ensures that the metric and curvature are well-defined everywhere on each Ej. We 
handle this excision by restricting the radial coordinate to the range r G (0, fmax]- Thus each spatial slice is taken 
to have an outer and an inner boundary. 

With these considerations in hand, we turn to a discussion of the functional differentiability of the action. This 
requires specifying what variations are to be fixed on the boundaries, and may require the addition of more boundary 
terms [21]. Variation of the action (34) with respect to H gives the boundary terms 


dr 


PnSil 


t — tf 

J t^ti 


lira 

e^O 


(37) 
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for some initial and final times ti and tf, and 


lim / dt 


02 


^|0| 


And variation with respect to Pq gives the boundary term 


lim / dt 
£-> 0 / 




We define the variational principle by fixing O at the end points by 

^{t, = a{t), n{t, e) = b(t). 


(38) 


(39) 


(40) 


With this choice the 50 variation is well defined. The surface term arising from the symplectic piece is zero because 
initial data and its subsequent evolution fix O at ti and tf. Lastly, to keep Pa free at the boundaries, we add a surface 
term to cancel the 5 Pq variation. The final gauge fixed action is 


S= dt 


dr (^PoO - Pphys) - lim f 


dt 


2 

M 


1 T — T max 


'MPn 


(41) 


The summary so far is that in the process of deriving this action, the dust field and its conjugate momentum have 
been eliminated from the system, and the remaining metric field and its conjugate momentum (O, Pn) describe the 
geometry. Furthermore, as noted in [15], the dust time gauge results in the conversion of the former hamiltonian 
constraint of pure gravity into a non-vanishing true hamiltonian. 

There are other possibilities for fixing the variational principle. For example, we could have gone without adding 
the second boundary term and instead fixed the momentum Po on the boundaries. The choice above is the simplest 
since it requires boundary conditions for only, and still permits a large class of interesting solutions. 

With the variational principle well-defined, the equations of motions are 



(42) 

sgn(f^) {^Pn - • 

(43) 

Physical conditions 



Let us consider the spacetime metric and physical properties resulting from these gauge choices. In fixing a gauge 
for the field variables, we obtained conditions which fix the lapse and shift functions (33). The resulting metric is 

/ n2 p2 \ 

ds^ = — f 1- j dt^ -I- sgn(n)n^Pn drdt + V?dr‘^ + r'^dO'^. (44) 

The metric is non-degenerate so long as 11^ > 0, which implies that sgn(r2) is constant. 


1. Deficit angle 

Consider the ratio P of the circumference of a circle with radius rg > 0 divided by the proper radius 

/ y/m dO 


T := lim 


/r=ro-Se 


J 

In flat space this is 27r, but in general we may have P = 27r — a(7’o) where a{ro) is a deficit angle given by 

/ \ 


(45) 


a(ro) = 27r 


1 - 


ro 


|••ro+e 

lim / |n|dr 


(46) 


e -)-0 
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The limit of vanishing radius tq, after first taking the limit e —> 0, gives a = 27r — jofiToyfj' This implies that 

constant-time slices of the above metric generically describe a conical geometry near the origin with a time-dependent 
deficit angle. Note that in 3D gravity, a negative deficit angle corresponds to a point source with negative energy 
[2, 6]. In order to have a positive semi-definite energy at the origin, one would require that 

|O(f,0)|>l. (47) 


2. Energy density 

The stress-energy tensor may also be written in terms of the phase space variables. From the action we find that 

(48) 


T,.„ := -4=1^ = 


y/—g 27r ^ 

We see that m is the time-time component of the above. From the metric, or using (32) and (12) with the positive 
choice of sign to solve for to, we obtain the energy density 


SttTm = V\9\m = |D|Pn - , 


(49) 


where the equations of motion (53, 54) have been used to simplify the expression. A positive definite energy density 
requires that the right hand side be greater than or equal to zero. The spacetime Ricci scalar is: 


= TO = 


Pn 1 

D2 


(50) 


Since we have excised the point r = 0 from the spatial manifold, the curvature scalar is missing a delta function 
contribution at the origin when there is a conical singularity. 


3. Horizons 

Congruences of future directed outgoing and ingoing radial null geodesics are 

These satisfy = —2, and provide the null expansions. In our context these are physical phase space observables 
which are potentially useful in a quantum theory [13]. The outward null expansion of circles embedded in a spatial 
slice with unit outward normal = (0,1,0)/D is 

e 4,.)V„.„ = i(f-T). (52) 

Dynamical apparent horizons are obtained by solving Q(t,r) = 0 to give the horizon radius rh{t). This may have 
multiple solutions on a given time slice (see eg. [22, 23] for explicit examples). 


IV. CLASSICAL SOLUTIONS 

In this section we discuss classical solutions to our model. We find a large class of exact solutions and provide 
several examples. In particular we obtain a static solution for A < 0 which describes the BTZ black hole, and for 
A = 0 we find travelling wave solutions. 

To this point we have left the sign of D arbitrary. As noted above, we must require that sgn(D) is constant 
throughout the spacetime in order for the action to be well-defined. This implies that the solution space is split into 
sectors with sgn(D) = ±1. To keep the presentation simple, we assume D > 0 for the remainder of the article. The 
solution space for D < 0 is nearly identical with only trivial differences. 
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A. A / 0 


For the case of non-zero cosmological constant the equations of motion are: 


tl = 

—n', 

2 ’ 

(53) 

Pa = 

^Ph- 2Ar. 

2 ” 

(54) 


The second equation is similar to inviscid Burger’s equation, but with a source term coming from the cosmological 
constant; it contains only the momentum and can be solved independently. This is coupled to the first equation which 
resembles the advection equation but with a variable speed of propagation given by Pn/2. As we will see, any solution 
for the momentum then determines how initial data fl(0,r) evolves. 


1. General solution 

There is an auxilliary, flat spacetime with Lorentzian signature defined by the (t, r) plane. On the auxilliary 
spacetime the momentum equation (54) is a conservation equation daJi = 0 for the current: 


Ji = 


which has an associated conserved charge given by: 


—Fb, 


p2 


(55) 


qi = lim 
£->■0 


Pq dr. 


Considering the system as a whole there is another conserved current: 


J2 = 


Pa 


n\2Ar- I -^P/, - 2Ar 


Pa ( Pa 


where the conserved charge is: 


<72 = lim 

e ^-0 


O ( ^Pa - 2Ar ) dr. 


(56) 


(57) 


(58) 


It is well-known that conservative equations can be solved by the method of characteristics. To employ this method, 
we consider characteristic lines parameterized by s, described in terms of parametric equations for the coordinates 
r(s) and t(s). Differentiating with respect to s, we obtain: 


dt 


, dr 


ds 


Pa — PaT: — Pa~^ — —2Ar. 


' ds 


ds 


This is equivalent to the equation of motion if we have the following equations along each characteristic: 


m 


= 1 , 


dr 

ds 


Eal 
' 2 ’ 


ds 


Pn = -2Ar. 


These equations are solved by: 

t = s, r = rQ cosh VTis — 


2y/A 


sinh VEs, 


Pq = Pq cosh y/Xs — 2\/Aro sinh VXs, 


(59) 


(60) 


(61) 


where the initial values are ro = r{s = 0) and Pq = Pn(s = 0). Each characteristic is labelled by the ‘starting point’ 
ro, and initial data for the momentum is a function of the radial points on the initial slice Pa(t = 0.,r) = Po(ro). 

Given a solution for the momentum we can solve (53) for G using the characteristic method again. The characteristics 
for this equation are the same as those for the momentum equation of motion, but here we have: 


^n = o, 

as 


(62) 


so that n is constant along each characteristic. This implies that given some initial data 17(0, r), the configuration 
field simply flows along the characterstic lines defined by the momentum. 







2. Examples 


Static solutions are obtained by setting = Pq = 0 in (53, 54). This gives 


n = Cl, (63) 

Pn = ±2V'C2+Ar2. (64) 

The metric may be put in a more succinct form as follows. Rescale C 2 and the r and 6 coordinates by absorbing 
the constant Ci as r = Cir, 0 = 0/Ci and C 2 = C^C' 2 . With this rescaling the angular coordinate has a range 
0 < 0 < 2-^1 Cl so that a choice of Ci > 0 implies a deficit angle defined by as described in the preceeding section. 
The line element becomes 


= —fdt^ ± 2y^l — / df dt + df^ + r^dO^, 


(65) 


where f {f) = \ — C2 — Af^. 

This solution remains well-defined for any choice of C 2 > 0, which is required for Pn to be real at each point. For 
de Sitter spacetime (A > 0) there are no additional restrictions, but for the anti-de Sitter (AdS) case (A < 0) the 
radial coordinate has a limited extent in order to keep Pn non-imaginary: 


0 < r < 



( 66 ) 


Let us consider the AdS case further. The above line element is in fact a generalization of the BTZ black hole which 
allows for a deficit angle due to the choice of Ci. This can be seen by transforming to a new time coordinate 


t = t ± 



\/l - f{x) 

fix) 


dx, 


which puts the line element in the form 


ds^ = -fdP + f-^df^ + f^de^, 


(67) 


( 68 ) 


where we note again that the angular range is 0 < 9 < 27r/C'i. 


This spacetime has an event horizon at f = 



1)/|A|, and when Ci 


1 it is the BTZ spacetime in flat slice coordinates. 


B. A = 0 


With zero cosmological constant the equations of motion have a remarkable symmetric form 

ri = -Pfifi', 

Pq = 


(69) 

(70) 


The momentum equation of motion is now Burger’s equation with vanishing viscosity. There is a substantial volume 
of literature on the subject. Most notably, this equation gives shock waves when characteristics cross. 


1. General solution 

The equations of motion are again conservation equations daj°' = 0 in the auxiliary flat Lorentzian spacetime 
defined by the {t, r) plane. With A = 0 the currents are 


Ji = 


—Pn, 


p2 


J2 = 


-n 


p2 


p2 

n-9-pf 


(71) 
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and the corresponding conserved charges are 

P'^Tna.x ^ P'^Tn&x 

qi = lim / Pn dr, q 2 = lim - / ^PqPq dr. (72) 

The equation of motions can again be solved by the methods of characteristics; one need only put A = 0 in the 
equations from the last section. There are two classes of solutions: 1) Pq is constant and fl = h(r + -^t); 2) Pq is 
not constant and O = /i (Pn), for an arbitrary function h. This last fact is a remarkable consequence of the structure 
of the A = 0 equations. When Pn is constant the characteristics are guaranteed not to cross. 

We also note that for vanishing cosmological constant the parametric equations (61) can be inverted to yield the 
following: 


2r + Pnt - 2/(Pn) = 0. 


(73) 


Given any choice of function /(Pn) of the momentum, solutions to (70) are given by the roots to this equation. 


2. Examples 


Let us note three types of solutions. The first is a class of static solutions obtained by setting Pn = 0 in (69-70). 
This implies = /(r), a nowhere-vanishing, but otherwise arbitrary function. The resulting metric is 


= -dt^ + /{r^dr'^ + r'^dO'^. 


(74) 


2 /' 

The energy density is given by 2'K%tt = —pr, a-nd its sign is determined by the sign of f. The spacetime Ricci scalar 

/ 


f 

is and the r~^ factor indicates that solutions are generally singular at r = 0, except for the particular 

choice /(r) = ±(Ci — C' 2 r^)“^/^. Constant time slices are cones with deficit angle a = 27r(l — l//(0)), and there are 
no horizons in this spacetime. 

The second is a self-similar solution obtained by setting /(Pn) = 0 in (73) to obtain 


Pn = -- 


2r 


(75) 


It is immediate that this solves (70), and leads to the equation of motion 


n = —o'. 

t 


One solution of this is O = 1, which gives the metric 




2r 


ds‘‘ = — ( 1-^ dP - drdt -I- dr^ + r'^dO'^. 


t 


(76) 


(77) 


The Ricci scalar from (50) is = 2/t^ so there is a spacelike curvature singularity at t = 0. Looking at the 
condition (52) we find horizons at r = —t. Constant t slices are flat without any deficit angle due to the choice of 
0 = 1 . 

A third class of solutions is obtained by setting Pq = 2v for some constant u € R. The O equation reduces to the 
advection equation 


O = uO' , 


(78) 


which has the general solution 


O = h{r + vt) = h{u) (79) 

for arbitrary h and no restriction on u, where u = r + vt labels each (straight) characteristic. If we choose O = C, a 
constant, we have a flat metric with deficit angle a = 2 tt{1 — 1/C) as described in (46). 

For a non-constant O, the u > 0 and v < 0 solutions describe respectively radially ingoing and outgoing profiles. 
The u > 0 wave metric is 


~ [l ~ + 2vll/drdt -|- Id/dr'^ + r'^dO'^. 


(80) 
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We note the following features of these “wave” solutions. The Ricci scalar from (50) is 



(81) 


and there are dynamical horizons if il,(t,r) = Ijv. Thus horizons will be present if the maximum value of the wave 
profile exceeds 1 /u, and the minimum is less than 1/v. 


From the expression for energy density (49) we see that the energy flux is positive where fl' > 0 and negative where 
n' < 0. There is a conical singularity at the origin with deficit angle a = 27r(l — l/fl(t, 0)). As the energy flux reaches 
the origin, positive flux adds to the deficit angle (which represents the mass of the singularity), and negative flux 
reduces the deficit angle. 


V. QUANTUM THEORY 


In this section we describe a non-perturbative quantization of the A = 0 theory. The circularly symmetric sector 
of the model we are considering has one local physical degree of freedom (Q), and as we have shown in the last 
section, the classical theory can be solved by the method of characteristics with the Pq solution providing a local 
(t, r)-dependent speed for the fl equation. The full quantum theory of this sector is more challenging. Although there 


is a physical hamiltonian and no constraints, the hamiltonian is unconventional in the sense that there is no separation 
of pure kinetic and potential terms. We can however achieve a full quantization of the Pq = constant sector of the 
solution space discussed in the last section. 

As we noted, this sector of the solution space describes either purely ingoing or outgoing waves. A first challenge is 
that since Pq is a constant, the symplectic structure we have been using up to now is not available. This is overcome 
by noting that we can obtain a new symplectic structure starting from the solution space of a differential equation 
[24]. The basic idea involves defining geometric structures on the solutions space that leads to a conserved symplectic 
current. The integral of this current on an initial value surface defines the desired symplectic form. 

For completeness, we now derive the symplectic two-form for the circularly symmetric, A = 0, Pq = 2v (where v 
is a constant), sector of the solution space and refer the reader to [24] for details. Once the symplectic structure is 
obtained, we move on to the canoncial quantum theory. 

In our case the differential equation is 


fl = vfl', 


(82) 


on the half plane t G (—oo, oo), r G (0, oo). For our purposes, we take this half-plane to define an auxilliary spacetime 
M with a flat Lorenzian metric given by 77 = diag(—I). 


Consider the space Z of solutions to (82). A point Vl G Z represents a solution to this equation, and a tangent 
vector do at this point is a small displacement which must also be in the solution space Z. Writing the displacement 
of the solution as Q -I- dQ, we find that the tangent vector 6Vl must also satisfy (82). 


The space of one-forms on Z is dual to the tangent space. If we label the spacetime points x G M, then the one-form 
dual to 5Vl is given by (5Q(a;). It is important to note that these one-forms are anti-commutative. 


5ri{x) 6Vl{y) -b SV,{y) 6V,{x) = 0. 


(83) 


The symplectic current is defined by 


Ja{x) = Sfl{x)daSfl{x) 


(84) 


which due to the equation of motion (82) and the anti-commutivity of one-forms, obeys the conservation equation 
y'^^daJb = 0. The associated conserved charge is given by integrating over a spatial hypersurface: 



(85) 


This conserved charge is the symplectic two-form we seek. It implies that the momentum conjugate to Q is 11 := O, 
with the equal-time Poisson algebra 



( 86 ) 


These are equivalent to the Poisson brackets for free scalar field theory, but the Hilbert space we define next will differ 
in that it includes only the ingoing or the outgoing modes, depending on the sign of v chosen. 
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Having obtained the symplectic structure, let us consider the Hilbert space we will use for quantization. Consider 
the positive energy (dust time) mode functions 


± 6*'=’') , fc > 0. 


(87) 


These sets have different boundary conditions at r = 0: = 0) = 2e and ip {t,r = 0) = 0. Both sets are 

orthogonal and complete on the half-line: 



dr ip^{r,t) ipk'{r,t) 



dk V;^(r,t) ipkir'.t) 


2TrS{k — k')^ 


27r(5(r — r’)^ 


and also satisfy 


( 88 ) 



dr ipk{r,t) ip^,ir,t) = 0. 


(89) 


Let Hy denote the Hilbert spaces with the bases ip ^, and let Hv = Hy © . The purely ingoing (outgoing) wave 

solutions are obtained by the normalized linear combination 


9 k{r,t) 


-^{iP+{r,t)+iP^{r,t)) 



^—ik(vt+r) 


(90) 


Clearly € Hv are solutions of our model. They may be viewed as “quasi-particles” if the bases given above for 
H^ and H~ are viewed as “particles.” The Hilbert space of the entire wave sector (all ingoing and outgoing modes 
labelled by u € M) is the tensor product 


H = ®vHv (91) 

Let us demonstrate the quantization in the v = 1 component. H(r’, t) and its conjugate n(r, t) may be represented 
as operators in Hi, in the manner that is standard in field theory: 


n(r,t) = 

/ {^kgk{r, t) + dlgk{r, t)^ , 

(92) 

n(r,<) = 

-i j dk Vk [akgk{r,t) - d\gk{r,t)^ . 

(93) 


Their commutator algebra implied by the Poisson algebra (86) leads to the usual commutators for ladder operators 

[a\,d\,] = [dk,dk']=Q [dk,d\,] = 5{k - k'). (94) 

Fock states are constructed by starting with the vacuum state |0fe) defined by afe|0fc) = 0, and the n^—particle states 
by 


\nk) 


1 



|0fc). 


(95) 


The Fock basis is given by products of the rife-particle states with different k values. 

This completes the specification of the quantum theory for the Pq = constant sector of the solution space. What it 
shows is that this non-perturbative sector of 2+1 gravity coupled to pressureless dust in spherical symmetry is exactly 
dual to the quantum theory of a massless scalar field on the half line. 


A. Metric operator 

With the above quantization we can now proceed to describe the “quantum geometries” for this spherically sym¬ 
metric sector of the model. The metric contains only the function H so it is possible to define the metric operator in 
the Hilbert space H by gab ■= ffob(f^)- A notion of geometry is given by the expectation value of this operator in a 
quantum state. There are thus an infinite set of possible geometries depending on the choice of state. 
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The metric contains the factor Ct'^ so we need to select an operator ordering of a^, to define it. This is provided by 
imposing the physical requirement that expectation values in semiclassical states give recognizable classical solutions. 
One choice for such states are the coherent states defined by 

ak\ak) = ak\ak) (96) 

These states are explicitly given by [25] 


\ak) = e 

nfc=0 



\nk)- 


(97) 


It is known that the expectation values of normal ordered operator in these states give the corresponding classical 
results. We therefore define 

gab ■= gab{- :) • (98) 


As an explicit example let us consider the state which is the vacuum for all modes except fc, and the coherent state 
for mode fc, 


i^)=i«fc)nio^-)- 

i/fc 


(99) 


This gives 

(: = l {igk{r,t)r al + (^^(r,f))^^^)^ + 2|a,H , (100) 

where ak are any complex numbers specifying a classical solution. The complex number ak must be such that the 
right hand side is positive definite in order to avoid a degenerate metric, and depending upon the value of ak there 
may be apparent horizons. The quantum fluctuations in these states 

A(fl2) = ( 52 . . ^ 52 .)2 


is of course not zero, since \ak) are the minimum uncertainty states. 

The expectation value of the metric in the rifc-particle state jrifc) on the other hand gives the metric 


ds^ = -(l- dt^ + ^drdt + + r^de\ 


k ) 


2nfc 


The constant time slices are cones with deficit angle 

a = 27r I 1 — 


k 

‘2Tlk 


( 102 ) 


(103) 


Recall that in 3D gravity, a conical singularity corresponds to a point source with a mass proportional to the deficit 
angle [2, 6]. This implies a discrete mass spectrum of the rife-particle states determined by the wave number k. It 
asymptotes to 27r as Uk gets large, and has a positive or zero energy for 2nfe > k. 

For special values of the parameters satisfying 2nk = k, the apparent horizon function vanishes everywhere 0(r’, t) = 
0. With these values, the conical defect of the spacetime slicing is such that the outward going null geodesics remain 
at constant radius. 

The self-adjoint metric operator defined above using the creation and annihilation operators provides, via the 
expectation value, a correspondence between quantum states and spacetime geometries. The coherent states lead to 
classical geometries with fluctuations. There is also the interesting possibility of obtaining “macroscopically entangled 
geometries” by using states that describe entangled superpositions. Construction of such states requires either two 
systems described by a tensor product Hilbert space (such as that for two spin one half particles), or the division of 
a Hilbert space into two subsystems. 

In the quantum theory of the wave sector we have discussed, it is possible to produce entangled states by considering 
for example states in the subspace of the full Hilbert space (91). Let |afc)i denote a semiclassical state in 'Hi 

and |/?fe)-i a semi-classical state in H-i. Then an example of a macroscopic entangled state of spacetime geometries 
is 

IV”) = (|afe)i|/3fe)_i -b |/3fe)i|afe)_i). 


(104) 
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There are numerous examples of this type involving two or more subsystems, even within a fixed v sector of the 
Hilbert space, but with states labelled by different k values. In a full quantum theory of gravity it would presumably 
be possible to divide the physical Hilbert space into sectors corresponding to, for example, black hole and cosmological 
geometries. It would then be possible to construct interpretationally challenging macroscopically entangled states. 


B. Quantum horizons 

From the forgoing we can consider the idea of a “quantum horizon” defined by a horizon operator [13] 

( 105 ) 

yZ 

which is the operator analog of the classical apparent horizon condition. In our gauge fixed setting this is physical 
observable. It is clear that the fluctuation of this operator is non-zero in a coherent state, and so the corresponding 
dynamical horizons are not sharply defined as they are in the classical theory. For a coherent state with an ak such 
that the horizon condition (h) = 0 is marginally satisfied, fluctuations of the metric operator can lead to uncertainty 
in whether or not horizons are present at all. 


VI. SUMMARY AND DISCUSSION 

We studied a new model for gravity in 2-1-1 dimensions. Unlike most of the existing literature on 3D gravity, the 
model has a local degree of freedom which manifests itself as a metric function in the dust time gauge. The resulting 
theory has novel aspects in circular symmetry. The equations of motion are simple yielding several interesting classes 
of solutions, including waves; the latter provide a “midi-superspace” sector of the solution space which is amenable 
to Fock quantization. 

The quantization provides some interesting and precise results. Among these is the observation that horizons 
fluctuate, which we showed using semiclassical states. It is natural to expect that this result goes over to four 
dimensions where it could inform issues such as the so called information loss problem in black hole evaporation. In 
particular, metric fluctuations imply that the separation of the Hilbert space into states which are strictly inside / 
outside the horizon, as is common in computing entanglement entropy, is an ambiguous procedure. Metric fluctuations 
further imply that the time of apparent horizon formation may be ambiguous for any choice of time coordinate. 

Metric fluctuations also inform the “firewall” issue, which at its core requires exactly null non-fluctuating horizons as 
a fundamental assumption. If a horizon is leaky because it has fluctuations, then it is clear that the central argument 
based on the impossibility of simultaneous perfect correlation between modes across a horizon on the one hand, and 
perfect correlation between early and late time Hawking radiation on the other, ceases to be an issue: no perfect 
null horizon, no monogamy problem. It is possible that horizon fluctuations are small for large black holes if the 
appropriate semiclassical state is sharply peaked. But if a firewall were to form, its accompanying back reaction on 
the metric would obviously lead to high curvature fluctuations, and in turn to large horizon fluctuations in the early 
stages of its formation. 

As a last comment, our quantization also demonstrates an exact duality between (a sector) of the model and 
1-1-1 free scalar field theory on the half line. This in turn is dual to fermionic theory via the well-known Bose- 
Fermi correspondence in two spacetime dimensions. This means that the quantization we have presented likely has a 
fermionic description. 

This work is a first exploration of the the use of dust time to study quantization of gravity. Natural extensions 
of our approach are to the 3-1-1 theory Bianchi models, spherically symmetry, and other reduced sectors, such as the 
Gowdy models. 
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